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ABSTRACT 

We consider the contribution of radially excited 2 3 Pi^ states to ip' production at the 
Tevatron energy. Production of these states via the conventional gluon fusion mecha- 
nism and via gluon and charm quark fragmentation processes is considered. We find 
that it is possible to account for the data on if)' production from the CDF experiment, 
by taking into account the decays of these 2 3 P 12 states into if>'. 
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The CDF measurement |]J of <//?/> and i/j' production at large transverse momentum 
(pr)i n 1-8 TeV pp collisions at the Tevatron is one of the most challenging experimental 
results of recent times. There is a huge discrepancy between the experimental data 
and theoretical predictions of conventional models, much larger than the uncertainties 
that are typical of such theoretical predictions. 

The production of quarkonia at large px is usually described by the parton-fusion 
mechanism in the framework of the colour-singlet model || [| |J. In this model, the 
wave-function at the origin for a S'-wave (or its derivative for a P-wave) quarkonium 
is convoluted with the cross-section for producing a heavy quark pair with the proper 
spin, parity and charge-conjugation assignments. There is an uncertainty in the de- 
termination of the wave-function parameters, and taken in quadrature with that due 
to choice of scale, an overall uncertainty of a factor ~ 3 is expected in the predicted 
normalisation of the quarkonium cross-section. Within this normalisation uncertainty, 
the colour- singlet model prediction agrees with the J/ip hadroproduction data from 
fixed-target and ISR experiments over a reasonable range of pt 0, H |§ ■ 

At higher energies where 6-quark production becomes important, the decay of b 
quarks is a mechanism that contributes to J/ip and ip' production, in addition to 
the parton-fusion mechanism of the colour-singlet model ||. In the CDF experiment, 
however, the secondary vertex information has been used to separate out the 6-quark 
decay contribution. The remainder are presumably produced directly by the fusion 
mechanism: therefore, the discrepancy between the large J/ip and ip' cross-sections 
measured by CDF and the theoretical predictions shows that the fusion mechanism 
grossly under-estimates the direct J/ip and ip' production cross-sections. 

The above discrepancy calls for a new mechanism of quarkonium production at 
large px- At high-energy colliders, where gluons and charm quarks are copiously pro- 
duced, the fragmentation of these particles into quarkonium states has been suggested 
as an important mechanism for large-pr quarkonium production [0]. Even though the 
fragmentation process is of higher order in the strong coupling constant, a s , it can 
dominate over the direct quarkonium production via fusion at large pt, where terms 
of the order of p T /mQ can easily compensate for the suppression due to the extra 
power of a s . The validity of this suggestion is borne out by the explicit computation 
of the fusion and fragmentation contributions for the Tevatron energy, presented in 
Refs. [§, [|, [L0|. In these papers, it was shown that the fragmentation contribution 
is the dominant contribution at the Tevatron energy; the largest contribution comes 
from gluons fragmenting into x's, which subsequently decay into J/ip's. Together with 
the fusion contribution it can explain the large cross-section measured by the CDF 
experiment to within a factor of 2-3, which is also the typical uncertainty in the pre- 
dicted normalisation. In Ref. ||, the energy dependence of the fragmentation and 
fusion contributions was studied as well. It was shown that at ISR energies the fusion 
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contribution dominates over the pt range of interest, although the inclusion of the 
fragmentation contribution helps to improve the quantitative agreement with the ISR 
data. 

In Refs. ||, H , the production of if)' by fragmentation was also studied. Unlike J/ if) 
production, if)' production does not involve the P-state contribution. Consequently, 
the fragmentation contribution is comparatively smaller. The total predicted if/ cross- 
section is more than an order of magnitude smaller than the cross-section measured 
by CDF. 

In this letter, we explore the possibility that radially excited P-state charmonia 



contribute via their electromagnetic decays to ip' production. Some authors [IT], [12 
have recently suggested the possibility of a large contribution of the 2 3 P 12 states to 
if)' production at the Tevatron energy. These 2P states have been predicted |TJ| to 
lie about 200-300 MeV above the open charm threshold, but their DD decay width 
is expected to be suppressed due to L = 2 phase-space as well as dynmical effects 
|TT| , fHfl . Consequently, the branching ratios due to the decays 2 3 Pi )2 — > 2 3 Si + 7 can 
be substantial, and, hence, the 2P states could potentially contribute a significant 
fraction of the measured if)' cross-section [] jTT], |12" 1. We present the first quantitative 
estimate of this contribution. 

The theoretical set-up which we need to compute the pt distributions has been 
discussed in detail in Ref. [§]. We will compute the contributions from the fusion and 
the fragmentation processes to the direct production of if)' as well as to that of the 2P 
states. The large-p^ production cross-section for the fusion process is given as 

^(AB ^ (2S,2P h2 )X) = 
dp T 

I dydx 1 x 1 G a/A (x 1 )x 2 G b / B (x2)- — — — v^j( ab ~^ { 2 S,2P lt2 )c), (1) 
J Zx\ — xt& at 

where the sum, in the above equation, runs over all the partons contributing to the 
subprocesses ab — > (25*, 2P 12 )c. G a /A and Gb/B are the distributions of the partons a 
and b in the hadrons A and B with momentum fractions X\ and x 2 , respectively, and 
x 2 is given as 

x x x T e-y - 2r 
2xi - x T ev y ' 

where r = M 2 /s, with M the mass of the resonance, s the centre-of-mass energy and 
y the rapidity at which the resonance is produced. 

2M T 2p 1 



1 9 AiviT ^PT 

x t = V x T + At = — 7 =-, x T = —/= (3) 

V s 



1 The contribution of the 1 £>2 state to the %/)' yield is expected to be small. For a discussion, see 
Ref. O. 
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The expressions for the subprocess cross-sections, dfr/dt, are explicitly given in Refs. || 
and [116 . 



The dominant fragmentation contribution comes from the fragmentation of gluons 
and charm quarks. It is computed by factorising the cross-section for the process AB — > 
(2S, 2P lj2 )X into a part containing the hard-scattering cross-section for producing a 
gluon or a charm quark and a part which specifies the fragmentation of the gluon (or 
the charm quark) into the required charmonium state, i.e. 

da(AB - (2S,2P lj2 )X) = £ £ dz da(AB - cI)Z)^ (2S , 2Pl i2 )(z, fj), (4) 

where c is the fragmenting parton (either a gluon or a charm quark) and the sum in 
the above equation runs over the contributing partons. D(z,fi) is the fragmentation 
function. The fragmentation function which is computed perturbatively at an initial 
scale of the order of the charm quark mass is evolved to the scale typical of the frag- 
menting parton which is of the order of Pt/z, using the Altarelli-Parisi equation. For 
the production cross-sections of the gluons and charm quarks we will use the lowest- 
order cross-sections, da/dt(ab — ► cd), just as in the case of the fusion contribution. 
The full set of initial fragmentation functions that we need to obtain the 2S and 2P 
contributions are : D g ^ 2S §, D g _> 2P [0, D C ^ 2S @ and D C ^ 2P H, ||. The final 



cross-section for the fragmentation process is given by a formula similar to Eq. |l] but 
with an extra integration over z, or equivalently over x 2 . 

As mentioned earlier, a major uncertainty in these calculations is due to the wave- 
functions. In the computation of the fusion contribution these appear in the subprocess 
cross-sections, whereas for the fragmentation contribution, the fragmentation functions 
at the initial scale are proportional to the wave-function factors. The wave function 
factor Rq for if)' is reasonably well-determined [g], but we have to make assumptions 



for the parameters that describe the 2P states. In Ref. |17j the fragmentation function 
for g — ► x is written in terms of two parameters Hi and H' 8 , where Hi is related to R[, 
and H' 8 is a parameter that describes the g — ► \ fragmentation via a 5-wave cc colour- 
octet state. We assume that the Hi for the 2P states is about the same magnitude 
as that of the IP states [I3|. The parameter H' s has a large range of uncertainty even 



for the IP states. On general grounds one expects the H'% parameter for the 2P state 
to lie in the same range Jll[j. Firstly the Hi and H' % contributions to any P-wave 



quarkonium production are expected to be comparable in size. Secondly the wave 
function factor Rq, which is the colour singlet analogue of H' s , is similar in size for the 
IS* and 2S states. Therefore we assume the same H$ parameter for the 2P state as for 
IP; whether or not they are exactly equal is immaterial in view of the large uncertainty 
in this parameter. The parameter values used are Pq = 0.49 GeV 3 ||, Hi = 15.0 MeV 
and H' 8 = 3.0 MeV JTJ |7|. We take the masses of the 2 3 Pi and the 2 3 P 2 states to 
be 3.95 GeV and 3.98 GeV, respectively In our computations, we use |2l[ the 
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updated MRSD-' parametrisations |22| for the parton densities in the nucleon, evolved 
to a scale Q 2 = /i 2 /4, where /i is chosen to be M T for the fusion process, and equal to 
Pt° = Pt/z for the fragmentation process. The fragmentation functions are evolved to 
the scale Pt/z. 

Finally, we come to the branching ratios for the radiative decay 2 3 P 12 — > 2 3 S\ + 7. 
The absolute value of these widths have been estimated by rescaling the known values 
for the bottomonium system |IJ. One gets r(2 3 P li2 -> 2 3 S' 1 +7) ~ 0.1-0.2 MeV. The 



corresponding total widths have been estimated, assuming phase-space and dynamical 



suppressions into open charm states, to be ~ 2 — 5 MeV each ||11|| . The resulting 
branching ratios for 2 3 P 12 — * 2 3 S' 1 + 7 lie in the range 2 - 10%. The results presented 
below are based on the optimistic value of the branching ratios = 10%. 

In Fig. 1 we present our results for Bda/dpT for ip' production in pp collisions at 
yfs = 1.8 TeV and integrated over a pseudo-rapidity interval \r]\ < 0.5, and compare 
our predictions to the data from the CDF experiment [ffj. B = 7.7 x 1CT 3 is the i/j' 
branching ratio into leptons. As in the case of J/i/j production ||, we find the g — > P- 
state fragmentation contribution to be dominant. The inclusion of the contribution 
of the 2P states increases the ip' cross-section significantly, and assuming a value of 
10% for the branching ratios this is capable of explaining the bulk of the observed 
cross-section. It should be noted of course that this choice of the 2P — > 2S branching 
ratios corresponds to the upper end of the estimated range. Within the lattitude of this 
range the cross-section can go down by a factor of 5. Further there is an uncertainty 
of a factor of 3 on either side in the production cross-section of the 2P states, arising 
largely from the Hi and H' % parameters. Consequently the predicted normalisation of 
the ip' cross-section can be anywhere between the experimentally observed range and an 
order of magnitude below. Thus it is premature to rule out the necessity of more exotic 
contributions, e.g. from the hybrid states suggested in fll|| . It is remarkable nonetheless 



that the inclusion of the 2P contributions can push up the predicted ip' cross-section 
by the required factor of ~ 30 within the lattitude of the model parameters. 

In conclusion, we find the 2 3 Pi i2 contribution to ip' production at CDF to be very 
important. We have computed the fusion and fragmentation contributions to these 2P- 
states and to the 25* ip' state, and find that with an optimistic choice of the 2P — ■> 2S 
branching, it is possible to explain the anomalously large cross-section for ip' production 
reported by at the CDF experiment. 



We are grateful to Frank Close for communicating his work of Ref. |TT[] to us before 
its public distribution. 
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Figure captions 



Fig. 1 The cross-section Bda/dpT (integrated over the pseudorapidity range —0.5 < rj < 
0.5) for the process pp — > ip'X as a function of pt at ^/s = 1.8 TeV. The data are 
taken from Ref. |l|]. The different curves correspond to the direct production via 
fusion (dashed line), the gluon fragmentation contribution (dashed-dotted line), 
the charm quark fragmentation term (dotted line) and the sum of all contribu- 
tions (solid line). 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409232vl 



